arXiv: 1506.08419v2 [cond-mat.quant-gas] 9 Oct 2015 
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According to the mean-field theory a condensed Bose-Bose mixture collapses when the interspecies 
attraction becomes stronger than the geometrical average of the intraspecies repulsions, gi 2 > 
511322 - We show that instead of collapsing such a mixture gets into a dilute liquid-like droplet 
state stabilized by quantum fluctuations thus providing a direct manifestation of beyond mean-field 
effects. We study various properties of the droplet and find, in particular, that in a wide range of 
parameters its excitation spectrum lies entirely above the particle emission threshold. The droplet 
thus automatically evaporates itself to zero temperature, the property potentially interesting by 
itself and from the viewpoint of sympathetic cooling of other systems. 
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The mean-field and the first beyond mean-field contri¬ 
bution, the famous Lee-Huang-Yang (LHY) correction, 
to the ground state en erg y of a homogeneous weakly- 
repulsive Bose gas read [l| 

E/V = {gn^/2){1 + 128^^/150? -fi ...), (1) 

where n is the density and a > 0 and g = Airh^a/m 
are, respectively, the scattering length and coupling con¬ 
stant characterizing the interparticle interaction. The 
LHY correction originates from the zero-point motion of 
the Bogoliubov excitations and is thus intrinsically quan¬ 
tum. It is also universal in the sense that it depends only 
on the two-body scattering length and not on other pa¬ 
rameters of the two-body or higher-order interactions. 
Quite naturally the experimental observation of this fun¬ 
damental beyond mean-field effect came from the field of 
ultra-cold gases , where the gas parameter na^ and, 
therefore, the relative contribution of the LHY term, can 
be enhanced by using Feshbach resonances Q . Note how¬ 
ever that the effect is perturbative; for na^ ~ 1 higher 
order terms and processes, in particular, three-body de¬ 
cay, come into play. A different situation is predicted 
for spinor gases where quantum fluctuations lift the de¬ 
generacy in the ground-state manifold @,0 or lead to 
quantum mass acquisition m- 

In this Letter we point out that in a Bose-Bose mix¬ 
ture the mean-field term and the LHY term depend on 
the inter- and intraspecies coupling constants in a differ¬ 
ent manner. Therefore, one can independently control 
them and make them comparable to each other without 
ever leaving the weakly-interacting regime. In partic¬ 
ular, an interesting situation, impossible in the single¬ 
component case, arises when the mean-field term, oc 
is negative and the LHY one, oc is positive. Be¬ 

cause of its steeper density scaling the quantum LHY 
repulsion neutralizes the mean-field attraction and stabi¬ 
lizes the system against collapse. The mixture can then 
exist as a droplet in equilibrium with vacuum without 
any external trapping jl2j| . This phenomenon naturally 
suggests a proof-of-principle experiment for observing the 
LHY quantum correction. The droplet can be prepared 


from currently available homo- and heteronuclear atomic 
mixtures by tuning the inter- and intraspecies scattering 
lengths into the unstable (from the mean-field viewpoint) 
region and by releasing the trap. We argue that several 
properties of the droplet are very unusual for ultracold 
gases and can have interesting implications. In particu¬ 
lar, we predict that the droplet can have a peculiar ex¬ 
citation spectrum containing only a continuum part and 
very few or no discrete modes (usually associated with 
collective excitations: surface waves, breathing mode, 
etc.) This means that starting from an ordinary finite- 
temperature trapped condensed mixture one arrives at 
a macroscopic zero-temperature object; excitations cor¬ 
responding to the continuum spectrum evaporate when 
the trap is switched off. The droplet can then be adi- 
abatically manipulated (for instance, be trapped again) 
and used as a bath for sympathetic cooling of other sys¬ 
tems. Another interesting property of the droplet is that 
its density can be strongly increased by going deeper into 
the unstable region and one can, in a controllable fash¬ 
ion, prepare an unprecedentally dense ultracold gas and 
study its optical properties, inelastic decay, etc. 

Let us choose the path-integral representation for 
the Bose-Bose mixture governed by the action S = 
f d^rdt£l^i(r, t), 4'J(r, t), 4'2(r, t), ^'^(r, t)] with the La- 
grangian density 
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where g,i is the chemical potential and ni(r) = |'I'i(r)p 
density of the i-th component, ga = Pnanjrai and 512 = 
2'Kay2,lmy. are, respectively, the intra- and interspecies 
coupling constants, m-r = mim2/{mi + m2), and we set 

n = i. 

In the mean-field approximation the system is stable 
when the quadratic form is positive defi¬ 

nite, which requires simultaneously 511 > 0, 522 > 0, 
and 5^2 < 911922- The classical ground state is char¬ 
acterized by uniform condensate densities rii satisfying 
= pi- In ths weakly-interacting regime quan¬ 
tum fluctuations are weak and one can expand C up to 




2 


quadratic terms in and '1'^*, where 4'' = 'I'i — ^/ni, 
thus arriving at the Gaussian path integral. The valid¬ 
ity condition for this approximation is simply na^ <C 1 
(here and after for qualitative estimates we omit sub¬ 
scripts assuming that the masses are of the same or¬ 
der of magnitude, mi ~ m 2 ~ m, and the same holds 
for the densities, ni ^ 712 ~ n, and scattering lengths. 
On ^ 022 ~ 1 012 ! ^ Cl). The equations of motion 
SC/6'9* = 0 and 6C/S9) = 0 give two Bogoliubov exci¬ 
tation branches [3, 13 | 
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where LOi = ^/gunik"^/mi + (fc^/2mi)^ are the Bogoli¬ 
ubov spectra for the individual components. 

The LHY correction is the zero-point energy corre¬ 
sponding to the Bogoliubov modes ([3]). The explicit ex¬ 
pression for its density reads 
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and can be rewritten as 






m 2 gi2 322^2 \ 
mi ’ giig22 ’ gii?^i ) ’ 

(5) 


where / > 0 is dimensionless. For m 2 = mi and for 
trivial cases, such as = 0 or 512 = 0, the two ex¬ 
citation branches (|3]) have the form of the usual single¬ 
component Bogoliubov spectra and the integral in Eq. o 
is analytic. In particular, for equal masses f{l,x,y) = 
X]±(l + 2/ i -I- /4:^/2. In any case, the 

main contribution to the integral in Eq. o comes from 
momenta of order 1 /^h = y/mgn. 

Let us now introduce Sg = gi 2 + pgiig 22 and dis¬ 
cuss the unstable regime when Sg is negative but small 
compared to gn > 0 and 522 > 0. The mechani¬ 
cal instability can be understood by diagonalizing the 
mean-field term J2i j=i 2 /‘^ ~ where 

A- ~ Sgy/giig 22 /{gii + 522 ) is negative a nd small, A 4 . Ri 
(gii + 522 ) 72 , n_ = (niy^ -I - n 2 v^) /Vgii + 522 , and 
n+ = [niy/^ - n 2 y/^)/Pgii + 522 - It is energeti¬ 
cally favorable to maximize n/_ and minimize i.e., in¬ 
crease both densities while preserving the ratio 712/711 = 
const = 511 / 522 - The instability also manifests itself in 

the fact that E^^- becomes complex for small momenta 
k ^ ^m\Sg\n. However, for k ^ 1/^h ^ \/m\Sg\n, 
i.e., in the region mostly contributing to the LHY term, 
both modes E^^- and Ek^+ are insensitive to small vari¬ 
ations of Sg and, in particular, to its sign. Note that the 
global increase in densities leads to the hardening of these 
modes and to the growth of the corresponding zero-point 
energy oc which is faster than the mean-field energy 


gain oc P. One can thus say that the long-wave-length 
instability is cured by quantum-mechanical fluctuations 
at shorter wave lengths. 

More formally, let us introduce a characteristic mo¬ 
mentum Pc such that m\Sg\n Pc -C 1/Ch- One can 
then obtain an effective low-energy theory by integrating 
out the modes with p > pc from the initial problem. The 
result is that the effective Lagrangian density (for the 
low-momentum part of 'I'i) is given by Eq. ([2]) minus the 
zero-point energy of the high-energy inodes. The latter 
is given by Eq. 0 or 0, where one can set 5^2 = 511522 
neglecting small finite-Jg corrections and extending the 
integration interval Pc<k<ootoO<k<oohy using 
the fact that the low-momentum contribution to Eq. o 
is negligible. 

As a result of the competition between the attractive 
mean-field term ex P and repulsive LHY term oc the 
mixture can exist at finite density without any trapping, 
i.e., in equilibrium with vacuum. Note that the LHY 
term competes with the attractive term A_n/. but it is 
still much too weak compared to A+n^. This locks the 

ratio of the equilibrium densities /n^i'^ = \/gii/g 22 . 
Then, omitting the numerical prefactors, the structure 
of the interaction part of the energy functional reads 
gn+ — \Sg\P_ + 777^/^((777_)®/^ leading to the equilibrium 

value 71 )°^ ~ 5g‘^/m?g^ ~ {5g/g)'^/P. This means that 
the system remains weakly interacting (gaseous param¬ 
eter np -C 1) under the condition [Sg/g)'^ -C 1. More 
quantitatively, we obtain 

„(o) = ^ ^ ^ (6) 

^ 1024 p{m2/mi, 1, y'g22/gii) «ii 5ii522 ’ 

and explicitly for equal masses 

( 0)1 _ 257r (ai 2 + 

P - 1024 aiia22»70IT + ^ ’ 

For finite particle numbers the system is in the droplet 
state. In order to study its density profile and low- 
lying excitations we simplify the problem by setting 

'I'i(r,t) = '\Jri'p(j)[v, t), where is a scalar wave 

function M- This assumption neglects possible relative 
motion of the components, i.e., energy-expensive fluctua¬ 
tions of 77_|_, and more exotic situations such as, for exam¬ 
ple, vortex excitations characterized by different charges 
in the two components. However, it is justified for the 
ground state of the droplet and the low-energy part of its 
spectrum. With this reservation, the effective theory for 
the field (j){r, t) is governed by the action 

S = r] J Prdi[Re{i(j)*dp) — + fl\cj)\'^], ( 8 ) 

where 5 = ( 2 / 3 ) 1 ^ 5177 ^°^ 77 ^°^ and e{(j),(j)*) = 

|Vf(/P/2 — 3|(/|"^/2 -I- |(()|^ is the rescaled energy density. 
We have introduced the rescaled coordinate f = r/^ and 
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time t = t/r, where 


/3 y^lmi + _ 3 

(9) 

In what follows we denote rescaled quantities by tilde. 

Since the coefficient r\ is large, t] ~ ^ 1, 

the problem dS]) can be treated (quasi)classically. The 
equation of motion reads 

= (-V?/2 - 3|<^p + 5I0IV2 - m (10) 

and the Gross-Pitaevskii (GP) equation for the ground 
state is obtained by setting (j){r,i) = (j)o{r) El- The 
chemical potential fl is fixed by the normalization condi¬ 
tion N = J (fir \(j)Q\'^, where N is related to the number 

of particles of the i-th component by Ni = The 

uniform solution for infinite N corresponds to c/jq = 1 and 
/i = -l/2. 

For large but finite N the ground state is a spherical 
droplet of large radius R « (37V/47r)^'^^ with approxi¬ 
mately unit bulk (saturation) density. Near the surface, 
if we denote by x the coordinate normal to it, the wave 
function (j)o{x) satisfies the one-dimensional GP equa¬ 
tion of the form (fi(j)o{x)/dx^ = —dU{(j)o)/d(j)o, which de¬ 
scribes the classical motion of a particle with coordinate 
(f>o and time x in the potential U{(f>o) [l^. The integral 
of motion {d(l)o jdx)'^ j2 + U{(j)o) = 0 can be integrated 
again resulting in the implicit formula for (f>Q(x)^ 
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from which we see that the surface thickness is of order 
^ (in initial units). That the corresponding momentum 
scale 1/^ is much smaller than pc a posteriori justifies 
the applicability of the low-energy effective theory ([8|). 
The surface tension equals a = J dx[e{(po, ((> 0 ) ~ d‘\4’o\'^] = 
3(1-1-V3)/35 and we obtain the spectrum of the droplet’s 
surface modes (ripplons) valid in the limit of large N [I^ 


= V^(4^/3)l(/-l)(l + 2)d/fV, (12) 

where ^ > 0 is the angular momentum. The dipolar mode 
(1 = 1) can be regarded as a center-of-mass displacement 
of the droplet, therefore wi = 0. 

In general, when N is not large, the size of the droplet 
is comparable to the surface thickness. In this case, we 
calculate the wave function (po (see the upper panel in 
Fig. 1) and determine the chemical potential fl of the 
droplet and its energy E = J d^re{(j)o,(j)Q) numerically. 
We also study its small-amplitude excitations by solv¬ 
ing the Bogoliubov-de Gennes equations obtained by lin¬ 
earizing Eq. m with respect to small ^(f,t) — ^o(f). 
We find that with decreasing N the droplet becomes 
metastable {E > 0) for N < 22.55 and unstable for 
iV < iVc ~ 18.65. This effect can be understood by 




Figure 1: Upper panel: the droplet wave function vs radial 
coordinate for N = Nc ~ 18.65 (solid), iV = 30 (dashed), 
N = 500 (dash dotted), and N = 3000 (dotted). Lower panel: 
the rescaled energy per particle E/N (dash-dotted), particle 
emission threshold —fi (thick dotted), monopole mode fre¬ 
quency Wo (solid), frequencies of higher angular momentum 
modes ui (dashed), and the corresponding surface-mode ap¬ 
proximation Eq. (II 2 II (thin dotted) versus (N — 


the following qualitative arguments. For a droplet of size 
R the energy is composed of the kinetic term oc N/R^, 
two-body interaction oc —N'^/R^, and LHY contribution 
cx With decreasing N the minimum of E{R) 

first becomes metastable (the global ground state cor¬ 
responds to i? —>■ 00 ) and then disappears completely. 
We should mention the analogy of this situation with the 
problem of a harmonically trapped scalar Bose-Einstein 
condensate with attractive two-body interactions (see 
for review). In that case the LHY correction is neg¬ 
ligible and the stabilizing role is played by the potential 
energy oc NR?'. There is also a metastable minimum 
which becomes unstable (in this case with increasing N) 
as a result of the interplay of the kinetic, interaction, and 
potential energies. 

In the lower panel of Eig. I we show the quantity E/N 
(dash-dotted line), particle emission threshold —fl (thick 
dotted line) which separates the discrete and continuum 
parts of the spectrum, frequency of the monopole mode 
Wo (solid line), and frequencies of higher angular momen¬ 
tum modes w/ (dashed lines) as functions of {N — Nc)^fi. 
The thin dotted lines extending above the particle emis¬ 
sion threshold represent the result of Eq. (fT^ . All ex- 
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citation modes cross the threshold for sufficiently small 
N. Only the monopole mode reenters at iV Ri 20.1. Re¬ 
markably, in the interval 20.1 < N < 94.2 there are no 
modes below —/r and, therefore, exciting the droplet is 
equivalent to spilling of particles or, more generally, to 
breaking the droplet into smaller pieces. We thus deal 
with an automatically evaporating object which is also 
macroscopic since the actual particle numbers are large, 
Ni = nf- {g/\Sg\f/^N » N. 

Let us now discuss excitations involving relative mo¬ 
tion of the components with finite n+, neglected in de¬ 
riving Eq. dS]). In the homogeneous case these modes 
correspond to E+^k and, for the droplet of size ^ 1/^, 
we can estimate the lowest frequency to be In 

rescaled units this quantity is approximately ~ 

\/ 5 / 1 ^ 5 ! ^ 1, i.e., all such modes are in the contin¬ 
uum. It is also possible to have a nonzero with¬ 
out exciting the relative motion. Indeed, starting from 
the configuration N 2 /N 1 = \Jgwjg 2 i and adding, say, 
excess particles to the first component makes N 
proportionally larger leading to a linear gain in energy 
{r]fT)fLSN ^ —\6g\nSNi. However, due to the quadratic 
penalty coming from the term A+n^i?^ ^ gSNf/R^ (here 
R = is the droplet radius) one can increase SNi/Ni 
only up to a critical value ^ \Sg\/g I. Beyond this 
point the excess particles no longer bind to the droplet. 

In the vicinity of the instability point the droplet can 
decay by quantum fluctuations similarly to the quantum 
tunneling of an attractive scalar trapped Bose gas to¬ 
wards collapse 0E1. Although in our case we are deal¬ 
ing with an expansion to infinity rather than collapse, in 
both cases the mechanism is associated with a soften¬ 
ing of the breathing (monopole) mode and its tunneling 
under a barrier. In order to estimate the corresponding 
lifetime we first note that near the instability point the 

chemical potential behaves as jl ~ jlc — y^2(iV — Nc)/N” 
[^ . where we find numerically fic ~ 0.061 and N" = 
Ri 2190. Let us parameterize (po by the 
chemical potential and introduce the derivative Sn{r) = 
\(j)o{f)\d\(j)o{f)\/dfL taken exactly at the critical point. An 
excitation “along” 5n{f) does not change N and corre¬ 
sponds to the monopole mode, the frequency of which 
vanishes at the instability point. Close to this point 
the tunneling path goes along this soft degree of free¬ 
dom and we introduce the corresponding coordinate A 
by writing the density as h{r,i) = \(j)o{r)\‘^ + X{i)Sh{r) 
and the wave function as (j){f,i) = sjn{f, i) exp[f0(f, i)]. 
Here the phase 9 is determined by minimizing ([8]) for a 
given trajectory A(t), which is equivalent to solving the 
continuity equation —n = Vf;(nVf0) with respect to 9. 
In this manner we obtain the effectively one-dimensional 
problem with the action 

S = vJ dt[meffA^/2 - [/eff(A)], (13) 

where the effective mass rUeS = 


TT J^[J^ r'^Sn{f')df']'^/[f^nc{f)]dr ks 2240 [here ndr) 
is the density profile at the instability point] and 
effective potential Ues{X) = (iV"/8)[(/lc — /i)A^ — A^/6]. 
The tunneling rate out of such a quadratic-plus- 
cubic potential in the quasiclassical regime equals 
r = Aexp(—5'®) [2l[, where the bounce exponent = 

a; 3.38g{N-Nc)^/'^, 
At = ujoy/lbS^/2 'k, and the monopole frequency 
ujQ = [N"{N — iVc)/8mgfj]^/"‘ (see Fig. [T|). Thus, for the 
observation of the macroscopic tunneling, due to the 
exponential dependence of T, the most relevant region 
of parameters is N — Nc ^ {^g/gY ■ 

Let us now discuss possibilities of having such droplets 
in experiments with ultracold gases. A very promis¬ 
ing candidate is the mixture of the second and third 
lowest hyperfine states, F = 1,771;’ = 0 (state I) and 
F = I,777f = —1 (state 2), of The corresponding 
scattering lengths Qij as functions of the magnetic field 
have been studied theoretically [H, and experimen¬ 
tally [ 13 . In particular, the condition ai 2 = —^011022 
is satisfied at Bq ss 56.77G where an ~ 84.3aBohr, 022 ~ 
33.5aBohr, and the quantity i 9 (ai 2 + v^aiia^)/9il ~ 
12.09aBohr/G. For example, for B — Bq = —250mG 
the saturation densities equal 77 ^°^ r; 3.3 x I0^^cm“^ 
and r; 4.9 x lO^'^cm”^, length scale ^ r; I.96/7m, 
and time scale r r; 2.4ms. Then, TV = 30 corresponds 
to TVi = 0.75 X 10®, TV 2 = 1.1 X 10®, and the droplet 
wave function is given by the dashed line in Fig. I (up¬ 
per panel). We can compare the characteristic time scale 
T with the lifetime rnfe ^ K^^{ni + 772 )“^ ^ I50ms, 
where we use the three-body recombination rate con¬ 
stant K 3 = 10“^®cm®/s [ 2 ^. It is also useful to note 
that T scales as whereas rufe oc \Sg\~‘^, i.e., there 

is still room for increasing \6g\ while keeping r < riife. 
Thus, one can have very dilute and slow droplets as well 
as dense and fast ones. Another interesting candidate 
is the dual-species mixture of ®^Rb and both in the 
F = I, mp = — I hyperfine state. This mixture is charac¬ 
terized by positive intraspecies interactions close to an in¬ 
terspecies Feshbach resonance at H = 117.56G, which has 
recently been explored experimentally [ 2 ^. In any case, 
unless we are dealing with the microgravity environment 
, the droplet has to be levitated against gravity. This 
can be achieved by using, for example, an optical poten¬ 
tial with a vertical gradient. Our theory can be modified 
to include a more general external confinement, say, har¬ 
monic. However, the self-binding and related properties 
make sense only if the system is unconfined at least in 
one (horizontal) direction. 
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